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ON ISOCHRONOUS DERIVATIVES OF THE FIRST AND 
SECOND ORDER IN SPACE DYNAMICS TASKS 


B. Ts. Bakhshiyan and A. A. Sukhanov 


INTRODUCTION _/3* 

In many problems of space dynamics, whose solution requires a linear ap- 
proximation, it is necessary to calculate the isochronous derivatives of the first 
order from the phase coordinates using their values at the initial moment of time, 
and also using the values of certain constants in the right sides of the equations 
of motion (such as the gravitational parameter, equatorial radius, etc.)* These 
derivatives are required when determining the trajectory of a spacecraft, for de- 
termining the reliability of predicting and the magnitude of the correcting pulse, 
for analyzing the perturbed motion and determining the erromof numerical methods 
of celestial mechanics [1-9]. 

The second isochronous derivatives may be necessary for a more effective 
solution of these problems, and also when considering the influence of nonlinearity. 

For example, they are necessary for studying the influence of nonlinearity upon the 
reliability of the least square method [10] and when analyzing the perturbed motion 
in the second approximation [11]. 

The matrix of the first isochronous derivatives satisfies a system of vari- 
ational equations and is related in a well-known way with the fundamental solution 
of the conjugate system [12], which in its turn plays an important role when opti- / 4 

mizing space maneuvers and flight trajectories [13]. 

In these problems, the presence of simple finite relationships are of great 
importance for isochronous derivatives and the solution of this conjugate system in 
the two-body problem. This is due to the fact that for slightlyj perturbed motion 
in several cases these relationships may be used to solve the problem with the 
necessary accuracy [2-8], and when it is necessary to consider the perturbations 

•k 
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they make it possible to greatly increase .the integration step of the variational 
system of equations (or the conjugate system) by using the method of Encke and 
Peano - Bekeiji^ 15 ]. it is desirable that these finite relationships have the 
following properties: 

a) convenience and simplicity for use on a computer, and great computa- 
tional speed; 

b) universality, i.e., the lack of singularities for parabolic and circu- 
lar orbits, and also in the case of rectilinear motion. 

The formulas for isochronous derivatives of the first and second order, 
given in [1,2, 4, 8 ] and [11] do not satisfy condition b) , In calculations using 
them for trajectories which are close to parabolic or "..rectilinear, large compu- 
tational errors arise. More complex formulas for the derivatives of first order, 
obtained in [3, 16] by differentiation of the universal solution of the two-body 
problem are suitable for any Kepler motion, but they do not yield relationships 
for the second derivatives which can be readily used on computers or the explicit 
expressions for the fundamental solution of the conjugate system. 

To calculate the isochronous derivatives, this article uses the relation- 
ship between the derivatives of the first order and the fundamental solution of 
the conjugate system (A. 1 and 2), In A, 3 relationships are obtained for deriva- 
tives of the second order. In A, 4 and 5 using the methods given in A, 1-3, formu- 
las are obtained for calculating the first and second isochronous derivatives:"of 
the two-body problem which are suitable for any type of orbits. Thus, for orbits 
which differ from parabolic orbits, the relationships obtained may be written in 
a form which is apparently the simplest for use on, a computer as compared with well- 
known methods, A . 6 gives formulas for calculating the first and second derivatives 
using the gravitational parameter of the two-body problem. 

A listing of all the formulas obtained is given in A, 7, 

We should also note that the integrals calculated’ in the Appendix are of 
independent value. 


The authors plan to publish a preprint in the near future, which will ob- 



tain the formulas given in this study in the form of a Fortran program. 


1. MATRICES OF FIRST AND SECOND DERIVATIVES 
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Let us examine a dynamic system described by the vector differential 


Cl.l) 


equation 

± = ^X, t) 

with the initial conditions (to) = 0 G 0 Here OC ^ 

*) 

is the. Jl -dimensional vector-column 

The matrix of the isochronous derivatives ^ (£) ~ • t^ 0 )°f the phase 

coordinates may be written in the following form, using their initial values 


<P 


2 x _ 


BXc L d x 


<5 X; 


,^ i • 


°y 


( 1 . 2 ) 


As is known [1], the matrix 
tions 


iK <p(t,to) satisfies a variational system of equa- 


_ <P =F<P , <P(t 0 ,tJ = I , 

wherejp ? , ; y I - 


( 1 . 3 ) 


unit matrix. 


The matrices of the second derivatives \ H =y H (t,t), K= 
may be determined as follows: " •" .-i 



' K 7 F 3 CCk , 

ijl'L dXddXo; J 


( 1 . 4 ) 


1 f • > ’ 

Below all the vectors will be assumed to be vectors-columns; the sign tifpi, 
designates transposition. ' 
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The study [9] gives the differential equations which are satisfied by the matrix 

K 


W' 


2. METHOD OF OBTAINING THE FIRST DERIVATIVES IN 
SOLVING THE CONJUGATE SYSTEM 

iy 

Let us look for the matrix i(^) as the solution of the system of equations 

A <£> - A© , (2.D 

where A = A(t) is a certain nondegenerate matrix of the ft.-th order; 
f ) 0 " A (to) . We shall use foy/W i to designate the rows of the matrix 
A , A 0 • Differentiating (2,1) with respect to time, taking into account (1.3), 
we find that \ f\^ satisfies the conjugate variational system of equations 

Ai ~ ~ A t F* , Ax (to) t Aoi ; i - i rr/ n . (2 - 2) 

Thus, the determination of the matrix is reduced to finding the fundamental 


n 


solution . r\ 1 ) of the conjugate variational system of equations. 


Let us assume that! {''^independent first integrals of the system (1,1) are 


known; 


/ — V 


Vt (x ) = V;(Xo) i = m . 

Differentiating (2.3) with respect to f'Y* at t = const, we obtain 

[ttO, 

3X 0 T ^ dXo 1 ) ? 


(2.3) 


(2.4) 




'0 w o 

It is thus apparent that the gradients of the first integrals of (1.1) 


A- - 


(2.5) 


are the ' HT row of the matrix A. In order to obtain the remaining ( ^ HT rows, 
it is necessary to findQ {resolutions of the system (2.2), which are independent 
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of each other and with the solutions (2.5)., For the case when the system (1.1) 
is reduced to the Hamiltonian form and n- 1 -. I of its independent first integrals 
are known, the last row of the matrix A may be found using the following theorem 
[12] (exclusion method). 

Charnyy method . Let us assume t t “T J. j * • •) Tl — i — gradients (2.5) 

of the independent first integrals of the system (1.1) is reduced to the Hamiltonian 
form. Then the n-dimensional row 


/8 


(Here 


to „ 

] 


( 2 . 6 ) 


0 1 

L-I o 

system (2.2), if the n-dimensional vector * 


is the matrix of nth order) is the solution of the 

and the function X/t) satisfy the 


following relationships : 


^ i jf ^ “ 1 ; • • • ; jX i ; 




(2.7) 

( 2 . 8 ) 


where C. are the arbitrary constants. Thus, the matrix A, whose row is 
1 • - I 

? l — I f • • • j ft , is degenerate then and only then, when the following 
condition is satisfied 

V: n Jl| 

S =£ 0 , (2.9) 

i=i 

I,-, 

where I yj . are the expansion coefficients of the Hamiltonian function gradient 
with respect to the rows 




3. OBTAINING THE SECOND DERIVATIVES 


Let us find the expressions for the matrices (1.4) when the fundamental 
solution of the system (2.2) is known in the form 

~ (x , x 0 ) t f to) t i — I,*** f ft. • 


/9 
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Let us write (2.1) using the elements of the corresponding matrices: 

!• • n •' . • 

pi ~ AoEl w; I = I f •• •} n • 


P=i 


Differentiating this with respect to Xoj , we obtain 

+ E P j ) <P fi +g Afp^ = D^+E^» 


(3.1) 


where 


F e - T) e .. F - 

X/p-7 , L-pJ , -Lfo LJ ) t—Olj. 


are elements of the matrix 


Txd SaI c BAT - rL 

D 7 >U> AM > 110 


•gj (3.2) 

8x ?.- dxo 9 y m: ’ 

A*! " 



Let us use QkC to designate the elements of the matrix 1 f\ . Multiplying 
(3.1) by KW , and summing the relationship obtained with respect to l£ , with 

r 1 ••-‘ni 

allowance for the fact that 

n 


a K e A ep ^/ = ^ 


e.P=i 


and changing to matrix form, we find the formulas for the second derivatives (1.4) 


^ K =ga Ki [DRE' -<f> T (E e +D l $)], K=i,...,n. « 


.3) 


We should note that if each of the matrices j £) / El " m ay substituted in (3.2) 

for certain \° Q in the fora of two components, so that the second components 

1 r i rr£ 


after substitution in (3.3) give an identical zero, then we shall use ^ ; E 
in this case to designate the first components. This situation occurs for the 
two-body problem (see A. 5). 


4. FIRST DERIVATIVES FOR THE TWO-BODY PROBLEM 


II Q 


In this section, formulas are obtained for the first derivatives which do 
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not have singularities for any types of Kepler orbits. 


For the two-body problem, the vectors 


xyf 


and the matrix F have 


the form 


X - 


where 


7 

V 


vM 


v • 

l -m 


F _fo I 

6 


(4.1) 


(4.2) 


where 7. V are the vectors of the coordinates and the velocities of the 

point, 7=|7| , V = /U| , ;U - — gravitational parameter of the 

attracting center. We first obtain the first five rows of matrix A, which are in- 
dependent of each other for any Kepler motion. We use the first integrals of the 
two-body problem for this: 


V C = 7 * 1 /, 

7 


r ; , IT* 


— n r 2 — 


K = v 


2 JU 


(4.3) 

(4.4) 

(4.5) 


Let us assume 


• 

are certain constant three-dimensional vectors. Multi- 


Pl > fz 

plying each of the equations (4.3) and (4.4) in a scalar manner by , > pi' 

and differentiating the relations obtained and (4,5) with respect to 7 ; Tv , 

according to (2.5), we find the first five rows of the matrix A: 

A; =[ 2jjj r , U-i 3 ^ _ _ j _ — 

A;, ? = [(\^* U ;+ V* Uf;) T , ( 


where 


U; = p ; x 7 , Xtrj =Vxpi, i = 1,2 . 


(4,7) 


For all types of orbits, except for rectilinear orbits, the rows ~ f\ 5 

are independent for any non-colinear vectors f P 2 * t * ie case 


/II 
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rectilinear orbits and orbits close to them, for the linear independence of these 


rows, it is sufficient,£for example, that the vectors f p2 be 

linearly independent of each other (this may readily be established by substituting 
U=a7 in (4.6), where Q is an arbitrary factor ), Therefore, for the uni- 
versality of formulas (4.6), it is advantageous to select jp jy. » , from this con- 
dition in every case. As Pi ^ p2 » we ma y select, : for example, two unit vectors 

of the system of coordinates corresponding to the components of the vector 
which are the smallest in terms of modulus. 


In order to find the row 


\^6 ■'» 


we may use the Charnyy theorem. We shall 


search for a vector 


'f 


which satisfies (2,7) and (2,8) in the form 


oh 7+jhU 
oC 2 T J 


(4.8) 


Substituting (4.8) into (2.8), with allowance for (1.1), (4.1) and combining the 



Solving system (4.10), we find 

d, = Q T 7 + 2£ a % 


(4.11) 
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d 2 +p,0 = Q T V , 

;A = 2 a T 7 + 3 e, +M, 


(A. 12) 


(4.13) 


(4.14) 


where ’.Of ; are arbitrary constants. We may readily determine that the co- 
efficients o(j2 an< ^ are £ nc l u ^ et ^ in (2.6) and (2.7) only in the form of the 

sum (4.13). This makes it possible to set in (4.8) 


Pi-.-O . 


( 4 . 15 ) 


Thus, the vector 


t- • 


found using the formulas (4. 8) , (4.11) - (4.15), satisfies 


the condition (2.8) of the Chamyy theorem. It may be readily seen that this vec- 
tor will also satisfy the condition (2.7) if the constants - C,. are related 

with the constants n C by the formulas 

• t ''l 

, C > = Pj- ’ 

= ( ox C) T (p; * cgy ; 

C5 = — £ i h - Q T -\|/ . (4.16) 

Since the row A,. is a -.gradient of the Hamiltonian function (i.e., in (2.9) 

s 0 at i^5 .As- = 2 ), then (2.9) for the nondegenerate nature of 

the matrix A in the case considered may be reduced to the inequality \ Cs^O • 
We set 


• r* 

a = £ 


(4.17) 


where *~2. is an arbitrary constant. According to (4.16), for the nondegenerate /13 
nature of the matrix A, the constants C C must satisfy the condition 


£ t h + £ 2 Y 2 ± 0 (\{/ = | Yl ) 


( 4 . 18 ) 


(this condition is satisfied for any Kepler motion, since h and i'/ljj. simultane- 
ously vanish)) ' 
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Substituting the values of 


found from (4.8), (4.11) - (4.15), 


r . * * 

(4.17) into (2.6), we find the sixth row of the matrix A: 


A 6 — Ag 4- £2 Ag 7 


(4.19) 

aJ =[-3d(i-t,)? r +v T , 

27 T -3(i-i)u T ] ; 

(4.20) 

/\^-[(ju7-20s)7 T +|i 


(4.21) 

use the notation 




7. -7 '■* V- - (7 .,7 


;i; 7; : 7 777+77 

:.: v + r~ + 

(4,22) 

. s = V T 7 =C 2 -yw7 

(C=|C|), 

(4.23) 

1 5 = y t [ ^ ~ 

1 . ? ’ 

-U)-ju[\dt : 

(4.24) 


It is shown in the Appendix that the integral in (4.19) is calculated in the form 
of a converging series : : •' 

,f " tyt- 77 , T 2 + V 3 + 


: 7c/t = 



</w ' 

■'V 


(4:; 2 5) 


" ' ‘ 




oo 


where 




2n+i 


26) 




Jr 


~ n/ '"7 flj 


T -• 


2_V 2 ■ ' 

7 JH 1 _ 

-t„)+ 7TU '" %TJ& 


and the control time, 


V/u 


7' = 


// 


7 = 


dr 

dT 

cf 2 ? 

dr 2 


7 T U ; ■ 

•% t 


= - 0(7 + 1 


(4.27) 

(4.28) 


(4,29) 


(4,30) 


714 
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According to .the Appendix, at d + 0 (the motion differs from parabolic motion), 
the integral (4.25) may be calculated using the final formula 



(4.31) 


where 


■ Z in 

(4.19). In the general case, for nondegenerate conditions on any orbits of::the 
matrix A, consisting of the rows (4.6), (4.19), according to (4.18), we could take 

E^Sianfh), £. z r= L . However, for a more effective calculation of the 

it is advantageous to select the constants from the following conditions: 


row A 


6’ 


at 


|hf > 


( £ >0 — the given small quantity) we assume 


= 0 


and at 


/K| < £ 


i.e. , 


if 


7 :° 2 i. 

the motion is parabolic or close to it, we set -0 Z = t In the 

first case of the row A^, we use the simple formula (4.20) to perform the calcula- 
tion; in the second case — we use formula (4.21). We should note that with this 
selection of the constants £ | # ' £ 2 formula (4.31) cannot be used. The 
series is calculated only for values of ^ which are small in terms of the modulus, 
which leads to their rapid convergence. 


5. SECOND DERIVATIVES FOR THE TWO-BODY PROBLEM 


/15 


In order to calculate the second isochronous derivatives for known matrices 
A and c£> , according to (3.3) it is sufficient to find the matrices * 

, determined from (3.2). We introduce the notation 

= Tpt , Pi = V r pi. r-A 

0 a Pez 

“Pez 0 ‘ Rgx 

L - Pis/ “pe* . 0 ' -J •-.•rv vf'. (5.1) 


Pi- 


. si 


mmK&m 


Qe = + vp 2 T - 2 p e v T , 


11 



Ri = "ip{ + Pe ^ 


Differentiating the rows (4.6) with respect to T . - TJ, we obtain — with 

•• j '* T~\ £ 

allowance for the notation in (4.2) and (5.1) — the first five ma trices’ : 


D e = 


D- 2 - 


D 5 = 


0 Pi 

.pro j 

Re +peG 



Vi' 


Ql 
-G o 
L 0 I J 


: 'Wq 

Re- 




• ? • -.y. 



i j 




ft , 2 


(5.2) 


Since the rows (4.6) do not depend on \T7 ") T , according to (3.2), 

to > yo, 

■ ° . ... 

E { = 0 , i = i , • • • , 5 . 

When finding the matrices D 6 , E , we shall assume that the selection of the 
constants £ j j 62 in (^ • 19) is performed using the algorithm used in A. 4. 
Differentiating (4.20), we obtain for orbits different from parabolic orbits 


(|hl> £#?.)’*: 


Df = 


-to 
21 


-3(f-Ul J 


6 w 
i 

Now let us find the matrices 


(5.3) 


e: = o 

D' 


then (| h U C/J/Vj 

CL is determined from (4,16) at 

■ 


for parabolic orbits and orbits close to 

We note_that, if the inequality Cg •=£ 0 

£ 1 = 0 — is satisfied for the given values of 

the initial phase coordinates \ r j Q 7 JV , then for a fixed constant |q , it is 

satisfied also in a certain vicinity of these values (in view of the continuity 
of the function (4.16) with respect to ' t 7 IT ). Therefore,vwhen finding 

• O f ^ O 1 • • • 

the derivatives of (4.21), we shall assume that the constant |Q does not depend 

f?o , Do , although it equals numerically (4,17). Let us write the 


on 


row 


(4.21), using (4.17) - (4.24) in the form 


/16 


Amv r ^ + 2 \) ¥ T R)7 T+ dj?'T;u t , (p T 7)r- 2 (y t r) v T J , 
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where 


t 

7=: r— - jj 

I 

y7rf 2R 

R - \ 7 01 

« . We first find 

Vi dOC, 


0 / c 


then (where Ve" 


designates the complete derivative, i.e., taking into account the dependence 

x(x ). As is shown in the Appendix, the integral R may be represented 
in the form 


where 


2 HR = \ p - \p o y 


(5.5) 


= i p( 7 , ir, t ) = 3 (t- to) y/^(7 T v)7 : + 2i 2 v, o . 

From (5.5) ; (£<, = ip j TJ 0 f t 0 ) V. Jr ;;. V- ' 


6) 


where 




T =£ T (i>v,t) = W T (^- 2 iyrR^) 


(5.7) 


(5.8) 


(This equation may be written as ^ 0 Differentiating (4.5), (5.6) 

with respect to 7 ,p and substituting in (5.8), we obtain 


j) = A 6 T + 3 h(t-t 0 )^ + 5 , 


(5,9) 


where A is the row (4.19), in which 


£ i--C 2 , S 2 = 2 , 


_[UxC] 

- L Cx^J 


(5.10) 


^ r _( 5 & + v 2 X7 T fj)7+(3/)7-27V+3^b] 

_ («7+'T j U j +(7 t U) 2 ) : Z-27 2 (7 t V)u (5 - n> 



/17 
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According to (2.1), 


0*$ - floe = 0 . 


Ig 1' Hoe — w * (5 , 12) 

It may be seen from (5,7) that at h » 0 the following equation must be satisfied 


/ \ 




(5,13) 




where *5; p., ( y ^ — 0 , i,e., according to A, 2, the row *3 p in this 

case must satisfy the conjugate variational system of equations (2,2), Substitut- 
ing •)) in (5,11) (according to 4,5), we obtain 



— *7 


(5.14) 


-T^(Tv)r- ( jirt- 3 (•7 t v) z )u 
( 3ju 7 +(Tvf)7. - 2 * Z *( t T u) U 

Substituting (5.14) in equation .'(212), where the matrix F is determined from (4.1). 

we see that actually satisfies the conjugate variational system of equations 

and, consequently, equation (5.13) is satisfied. It may be readily seen that 


f -(t t U r ) 7 -2 7 for To; 


5 = 'S, + h 

We find from (5.7), (5.9), (5.10), (5.12), (5.13), (5.15) 


7 2 7 


(5.15) 
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where 


i;dx 


0 / c 




' 3 ( t - to ) Ux C - ( 7 T ly7-2? V 


07 + 7 2 7 


i ‘j 


(5.16) 


(5.17) 


Ic = l ( 7 0 , V o , t 0 ) . ■ , _ ... ^ 

Differentiating (5.4) with respect to M k 7 0) ,Vo 

and assuming that the 

vector \ “X^f is constant, in agreement with the note made before the relationship 
(5.4), we obtain 

Df = D f E£, , (5.3-8) 


14 



where, with allowance for the notation in (4.2), (4.22) = (4.24), (5.17), 

• 6 r2sG +yu7 1 +vij7^:S;I - ?? T ] ; [ 

• ~ L , . si - t f A&A i 


D,V - 


\) 7 


■'***: l V f’J 

n, .1 




*7^-T-. i;l>. ./v: ' - A ' -A.K i " 
C r.-v’ V’jr . 1 .’ i LV--i : • •' 


■ T : .V 

>0 .* 


(5.19) 


(5.20) 


(5.21) 


We should note - that 'the marricesjj)^ -■ .\ ^ pT equal the derivative (3.2) within an 
accuracy of terms which, according to (5,12), (5.13), may be disregarded in (3.3) 
(see A. 3) . 


6. DERIVATIVES WITH RESPECT TO THE GRAVITATIONAL 
PARAMETER IN THE TWO-BODY PROBLEM 
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In the two-body problem for isochronous derivatives of the first order of 
the vector of- the phase coordinates . • -QC ' witlj respect to the 

gravitational parameter } the following formulas hold 

_ _J_ / = _ cfa T i (6-1) 

dju 3jU' J - ^ X »A / 

[where ! | (^!) is the matrix of the derivatives (1.2). 

Relationship (6.1) is a result of the integrals obtained in [7] for vari- 
ational systems of equations. The study [5] found these relationships by integra- 
tion of the equations for the desired derivatives. 


Differentiating (6.1) with respect to OCo anc * ’ we obtain the formu- 
las for the derivatives of second order 1 


d z x 


d<P _ 

djuA: : 


mm 

x:‘w 6 


i; 


( 6 . 2 ) 
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3 2 X 1 In 3x , d<$> ^ ) 

" " ~ 3ju ( 1 dju + dju X 7 


(6.3) 


where 


\T/ • J.'-. V; n\ are the matrices of the derivatives (.14) 

L , ( r» = . If •••; 6 ) 


7. SUMMARY OF FORMULAS 

Let us calculate the formulas which make it possible to solve the conjugate 
variational system of equations and to calculate the isochronous derivatives of the 
first and second order for the two-body problem. 


The matrix of the isochronous derivatives of the first order (1.2) may be 
found from the system of linear algebraic equations (2.1), where 

A = A , A, = A(f 07 Vc, t B ) 

are the matrices of 6th order. The first 5 rows of the. matrix A, with allowance 
for the notation in (4.2), (4.3), (4.7), are determined from (4.6). As the three- 

dimensional vectors JDj j Pi'*' 11 an d (4.7), we may select two unit vectors 

of the system of coordinates corresponding to the components of the vector 

which are the smallest in terms of modulus . The 6th row of the matrix at 

jhl>£/V7 0 (motion is different from parabolic motion) is determined from 

(4.20) , and atj 7 0 ^(parabolic motion and motion close to it) — from 

(4.21) . Here is the energy integral; e > o - a given small quantity. The 
quantities in (4.21) are determined from (4.22) - (4.30). 
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The matrices of the isochronous derivatives of second order (1.4) for the 
known matrices P ^ (^7 are found from (3,3), where are elements of the matrix 

f\ — . For : ^ = J. ; • ; 5 > the matrices j' J) ^ in (3,3), with allowance for the 

notation in (4.3), (5.1), are determined from (5.2), and = Q,. At 

|n| > ^/7 0 , where £ \> 0 , the same holds for the first derivatives, and 
the matrix is determined from (5.4), SQ , and at |hf 46 ju/l 0 
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the matrices j E — with allowance for (4.2), (4.4), (4.22) - (4.24), 

(5.18) — are determined from (5.19) - (5.22). The integral in (4.24), deter- 
mined from (4.25) - (4.30), can be assumed to be already calculated when the first 
derivatives are found. 


2 


The matrices T") ^ ^ . t ~~ 1 6 in (3.3) represent the 

X*/ 0 9 - 1 — Of. l . 

values of 


D f , E* 


at 


7 = 7 C , V=Vo, t = t, . 


(7.1) 


^ P /i;/7 t ^(parabolic motion and motion close to it) , the 


Thus, for the case 

sum + tL 0 in (3.3) equals the matrix (5.20) under the condition (7.1), 

since the sum of the matrices (5.21) and (5.22) equals zero under this condition 
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The first and second isochronous derivatives of 

spect to the gravitational parameter -JU for known matrices 
..., 6 are determined from (6.1) - (6.3). 


7 , 7 / 


with re- 

• *■ .!• 1 
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APPENDIX 
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To calculate the integrals 


: R =$Tdt 


(A. 1) 


R = J dt (A. 2) 

to - 

we change to a new variable — the control time , determined from the differen- 
tial equation [4] 


dT _ JE 

ST ~ T 


, TfU =0 . 


(A. 3) 


We may readily obtain formulas (4.29), (4.30) from (4.22) and (A. 3) for the first 

and second derivatives of .*" 7 : with respect to . We find the following from 
(4.30) 


n-i \ t 


7 f2n -°=(-o(r 7 

rr (2n) = /_ 


r? f K ) _ S T • i.v /. , 


' n>i 


(A. 4) 


where 

(4.27). Substituting 


..r-v-.n-: ■ 

— is determined by the formula 


r-Jt from (A. 3) in (4.30) and taking 

the integrals from both sides of the equation obtained, we find (4.28) for \T . 


From (A.l) and (A. 3), we have 


) = jL f 
' .;*j t 


Expanding the value R (0) = o of the function 

st the point jj “7 , we obtain 


(A. 5) 


in Taylor series 
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oo 


y -1 1^)1 7_rr\ r \_ 

jju.b; ni y l > r 


l 


where 


= ^(' z2 ‘ r -' zr ' r ^- f 


<~2 , 77''+ I'o-i 


3 


,2 - 

T 5 



co 


n=2 


. 2n+l. ■ U ( > >' J(2ri)! 

We find the derivatives in (A. 7) from (4.30) and (A. 4) 

2V^ 2n ‘^- Vl/'K _ 7 J ; 


(A. 6) 


(A. 7) 


•• K =0 






n-l 


n 0 2.1-2 

n z - /r-r//V— /‘•'K 


2n ■ l 


K=i 







(A. 8) 


' - 1 ( 2 ^ - 2 ^ 12 * 12 ' -^ 4 ); «.» 

Substituting (A. 8) and (A. 9) in the series (A. 7), we obtain the formula (4.26) 
for ^ Q . We may show that the series (4.26) converges absolutely for all \T . 
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For the case when jOtTr U (the motion differs from parabolic motion) 
we calculate the integral (A.l) in final form. We should note that when the equa- 
tions (4.23) and (4.25) are used for the focal parameter! |3 the following rela- 


tionship holds*. 




tv t»\ 2 . 


,-J? 


Z „,z 


= 21-dT-T . 


(A. 10) 


From (4.24) 


• IT. 


,1 . 7 . ■ ■ 


T 


j 7 a c/3r = 7 7 7o + o( J 7 2 c/f — J 7 cfc’ , (a.h, 


and from (A. 3) 
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Taking the integral from both sides of the equation (A, 10) considering A. 11) and 
(A. 12), according to (A. 5), we obtain the final expression (4.31) for the integral 
R. 


From 


We now find the integral (A. 2). From (A. 2) and (A. 3), we have 

R = 

and (A. 3), we obtain 




7=U 


s' = dl i 7 u 


Using (4.4), (4.3), (4.29), (A. 14), we may obtain the formula 


(A. 13) 


(A. 14) 


— 7 = 7- 

•7’7'-o/77 . 


J from (4.4), according 

to (A. 3), (A. 12), we 

have 


■ ■ • „ t • 


: [iv 

Wt* C - [ 7 

c/T" = 





?o)*C-§(t 

-t) 





| Using the dependence (4.20), we calculate the integral 



r T - 

i] 7 7 m 

-I"? * 

A . ‘ — 


(A. 15) 


(A. 16) 


(A. 17) 


Taking the integral of both sides of the equation (A. 15), considering (A. 12), j 
(A. 13), (A. 16), (A. 17), we obtain: 



Substituting (4.3) into (A. 18) (we must substitute 


G=' VU) 

/ 


I 


in the product ry X C- )»(4.27) and (4,29), we reduce (A. 18) to the form (5.6), 
L 0 ■ v 

(5.7). 
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